Avoiding dark states in open quantum systems by tailored initial correlations 
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We study the transport of excitations on a V-shaped network of three coupled two-level systems 
that are subjected to an environment that induces incoherent hopping between the nodes. Two of 
the nodes are coupled to a source while the third node is coupled to a drain. A common feature 
of these networks is the existence of a dark- state that blocks the transport to the drain. Here we 
propose a means to avoid this state by a suitable choice of initial correlations, induced by a source 
that is common to both coupled nodes. 
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I. INTRODUCTION 

The transport dynamics of energy or charge in many 
physical systems can be described by considering trans- 
port on a network of coupled 2-level systems. Depending 
on the physical system, the transport can range from 
being purely coherent to being purely incoherent. In the 
former case the dynamics can be described by continuous- 
time random walks 1 , while in the latter case the dynamics 
follows from Schrodinger's equation, which for complex 
systems and certain choices of the Hamiltonian defines 
so-called continuous-time quantum walks^i Coupling of 
a quantum system with purely coherent dynamics to a 
bath of harmonic oscillators (for instance phonons) can 
lead to a mixture of coherent transport and incoherent 
hopping induced by the environment. The dynamics (of 
the reduced density matrix) can be described by a quan- 
tum master equation of Lindblad type^, where a certain 
choice of Lindblad operators defines so-called quantum 
stochastic walks^. 

In this work we consider transport on a V-shaped 
trimer configuration, resembling for instance a system of 
coupled quantum dots^~— . The transport will be gener- 
ated by connecting a (incoherent) single source to the end 
nodes of the V-shaped trimer and a (incoherent) drain to 
the middle node, see below. Most theoretical works have 
focused on the case where the source creates an excita- 
tion on a single node of the network^^. One can also, 
however, consider the possibility where a source is con- 
nected to multiple nodes of the network. Such a source 
can then create an excitation that is in a superposition of 
these nodes, leading to initial correlations between these 
nodes. 

A common feature of these V-shaped, or circu- 
lar, trimer configurations is the existence of a dark 
statei^iii - — . Such a state causes the excitation to be- 
come trapped in the network and therefore leads to a 
blocking of the transport in the purely coherent case. To 
overcome this problem, one can either introduce an en- 
ergetic disorder on the nodes^iii or couple the system to 
a suitable environment where the decoherence process 
destroys the interference effects that lead to the dark 
stated. Here we propose a third method: a suitable 
choice of initial correlations that are induced by a single 



source which is coupled to the end nodes of the network 
creates an initial state which is orthogonal to the dark 
state. This then causes the absence of the dark state in 
the transport even in the purely coherent case, leading 
to a complete transport to the drain. 

The paper is organized as follows. In Sec. |TT] we in- 
troduce our model and provide a detailed discussion on 
the exact mathematical implementation of a source that 
induces initial correlations. In Sec. IIIII we discuss the 
transport efficiency of different initial configurations with 
the help of both analytical and numerical computations. 



II. MODELLING THE TRANSPORT 

1. Coherent and incoherent dynamics 

We consider the dynamics of excitations on a trimer 
network. The coherent (quantum) dynamics on this net- 
work is described by the Schrodinger equation, or equiv- 
alently by the Liouville - von Neumann equation, with 
the general Hamiltonian 
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Here, Ek is the site energy of node \k) and Vki are the 
transfer rates between nodes \k) and \l). Note that for 
certain choices of the site energies and the couplings there 
exists an eig; enstate \D) = (|1) - \2))/y/2 of H , having 
only an overlap with nodes 1 and 2£ & n i 12 . 

Now, if the system is in contact with an external en- 
vironment, the total Hamiltonian takes the form H to t = 
Ho + H# + Hi n t, where Ho is the Hamiltonian of the 
network, H# is the Hamiltonian of the environment and 
Hi n t specifies the interactions between the network and 
the environment. When the environmental correlation 
time is small compared to the relaxation time of the sys- 
tem, one can describe the dynamics on the network by a 
master equation in Lindblad form 4 : 
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FIG. 1. An illustration of the two different ways of connecting 
the source(s) to the trimer network: (I) corresponds to the 
configuration with 2 independent transitions to the network 
and (II) corresponds to the configuration with a source that 
creates initial correlations between nodes 1 and 2. 



with the constants Xki > for all k and / and 

V(L kl ,p N (t)) = L klPN (t)Ll - i [LlL kl , p N (t)} . (3) 

For our model we assume the Lindblad operators to 
be given by ~Lki = \k) The term in Eq. (|2]) corre- 
sponding to the operator Ljm models the incoherent ex- 
citation dynamics, induced by the environment, between 
the nodes I and k with rate Xki- 

In the purely incoherent limit we assume simple hop- 
ping dynamics to be described by a Pauli master equa- 
tion (for the diagonal elements of p(t) only). The cor- 
responding transition rates A^, for k ^ /, can be phe- 
nomenologically estimated with Fermi's golden rule, e.g. 
Xki = A|14/| 2 . Here A is a constant that captures the 
particular details of the environment (for simplicity, we 
assume that A = 1). We further assume additional pure 
dephasing, which is taken to be identical for all the nodes 
of our network, i.e. we have Lindblad operators with 
\kk — A for all k , where A is the global dephasing ratei 5 -. 

Since the dissipative terms in the Lindblad master 
equation induce incoherent hopping between the nodes, 
we can introduce a parameter a, with < a < 1, that al- 
lows us to interpolate between purely coherent dynamics 
(a = 0) and purely incoherent dynamics (a = 1): 

= (1 - a)C coh (p N (t)) + a£ env (p N (t)), (4) 

with Cc6h(pN(t)) = -i[Ho,pN(t)] and C e nv(pN(t)) = 
£incoh(piv(£)) + £deph(pN(t)). This approach is also 
known as the quantum stochastic walk 5 . The genera- 
tor Ancoh corresponds to the terms in (|2]) that generate 
incoherent transfer between the nodes and the genera- 
tor £deph corresponds to the terms that generate pure 
dephasing. 



2. Sources and drains 

A source is included in our system as an extra node |0) 
that is incoherently coupled to the nodes of our network 
in order to make sure that there is only transport from 
the source to the network and not back. In general, tran- 
sitions from the source to a general state = a& \k) 
of the network can be phenomenologically modelled by 
the Lindblad operator L s = (0|, leading to the follow- 
ing additional term to the master equation: 

^source (psN(t)) = rV(L a ,psN(t)), (5) 

with T representing the rate at which the excitation flows 
into the network. Note that for a A^-dimensional network 
with a source, the reduced density matrix psn is repre- 
sented by a (N + 1) x (N + 1) matrix. For an initial 
preparation in the source node, i.e. Psn(Q) — |0) (0|, one 

can show that the density matrix can be written in the 
form i6. 

^)=( A ; (,) ;,))' 

where pNif) is the density matrix corresponding to the 
network nodes. 

In a similar fashion, we include a drain by coupling the 
state \N + 1) to the network. The incoherent transition 
from a state of the network to the drain with rate 7 
can then be modelled by the term 

(pSND(t)) = ^V(L d , PSND(t)), (7) 

in the master equation, with Ld = \N -\- 1) Thus the 
final reduced density operator psND(t) = pit) is repre- 
sented by a (N + 2) x (N + 2) matrix. 



3. Creating initial correlations with a source 

There are now two interesting ways in which we can 
connect the source to the end nodes |1) and |2) of the 
trimer network: 

(I) The source can either feed node 1 with rate T/2 or 
node 2 with rate T/2. We assume these processes 
to be independent of each other. We can model this 
with two dissipators representing the two indepen- 
dent processes: 

4oLe(PW) = (0| ,p(t)) + ^V(\2) (0\, P (t)) 

= §Poo(*) (|1> <1| + |2> (2|-2|0) (0|).(8) 

(II) The source can also feed a superposition state 
between node 1 and node 2, which can in general be 
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FIG. 2. Dependence of the EST rjn(a) on the angle cj) for A = 1 and Y = 0.5. On the left we show the general result for all 
values of cj) and a. The right figure shows the cross-sectional curves for cj) — 0,7r/2 and <j) — 7r, together with the result for 
ill (a). 



written as |^) = (|1) + e i(f) |2)) /y/2. We can model 
this process with one dissipator: 



4ZMt)) = rv(W (o\, P (t)) 



:Poo(t) 



|1)(1| + |2) <2|-2|0> 



-e- i<#, |l>(2| + e i</, |2)(l| 
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See Fig. [T]for an illustration of these two configurations. 
The key difference between these two choices is that 

(2) 

£source(p(£)) creates initial correlations, depending on 
the phase (j>, between the two nodes, while £source(p(£)) 
does not. How these initial correlations effect the trans- 
port properties will be adressed in the following section. 



III. TRANSPORT WITH AND WITHOUT 
INITIAL CORRELATIONS 




FIG. 3. The dependence on the dephasing rate A, with F = 0.5 
and 7 = 1. The red curves correspond to an uncorrelated 
source and the black curves correspond to a correlated source 
with phase = 0. 



In previous work we used the expected survival time 
(EST) 77 as a measure for the transport properties of the 
excitation in the network^. Here we use it specifically 
to study the effects of initial correlations on the trans- 
port. This EST is defined as the average time it needs 
for the excitation to move completely from the source to 
the drain: 

oo 

77(a) = J dt (1 - p N+ i yN+ i(t,a)) . (11) 




The following representation of the EST in terms of the 
Laplace transforms pkk(s) of the components of the den- 
sity matrix, allows for a more convenient way to obtain 
analytical expressions for the EST 16 : 

N 

77(a) = limV^(s,a). (12) 

s— )►() z — ' 

We denote the EST corresponding to a source feeding 
independently the two nodes 1 and 2, Eq. (jSJ), as 77/ (a) 
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and the EST corresponding to the source feeding into an 
entangled state, Eq. (j9]), as 7777(a). To illustrate the key 
effects, we assume, for simplicity, that E\ = E2 — E3 — 
1, V13 = V23 = I, 7 = 1 and V12 = 0. That is, we focus 
on the situation when there is no bond between nodes 1 
and 2. For these parameters, it follows from Eq. (fT2j) 
that the ESTs 777(a) and 7777(a) take the form 

m (a) = l/r + f(a) /g(a) (13) 
mi {a) = l/r+ [J (a) - h(a) cos 4>}/g(a), (14) 



(15) 



with h(a) = 4(1 — a) 2 and 

/(a) = 4 + a(17 + 13A) + 2a 2 (A(A - 8) - 19) 

+3a 3 (ll + A(9 + 2A)) 
g(a) = 4a(2 + A) - a 2 (15 + 7A) + a 3 (ll + A(9 + 2A)). 

The dependence of 7777 (a) on the phase (j) is therefore pro- 
portional to coscj). Its amplitude —h(a)/g(a) is a mono- 
tonically decreasing function of a and vanishes when 
a = 1. Therefore 7777(a) converges to 777(a) when a — )• 1, 
where they both reach the value 4 + 1/r. This also shows 
that for (j) e [0,7r/2) and (j> e (37r/2,27r], 7777(a) < 777(a) 
and that the converse result holds for <fi G (tt/2, 37r/2). 
In the limit a — >• we find that 777(0) = 00, while 
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To illustrate these analytical results we show in Fig. 
[2j for A = 1, the dependence of the EST 7777(a) on the 
phase (j) and compare it to 777(a). One clearly observes 
the cos ^-dependence of 7777(a), see left panel, and an 
infinite EST for a = and <j> ^ 0. The infinite EST can 
be understood by noting that for these values of a and 
(j) the dark state is not influenced by the drain and is 
a stationary state of the system, causing both EST's to 
diverge when (f) ^ 0. For = the state is orthogonal 
to \D), causing the absence of the dark state in the full 
dynamics and leading to complete transfer to the drain. 

In Fig. [3] we show for <\> = the dependence on the de- 
phasing rate A. We observe that for increasing values of 
A the EST 777(a) decreases, leading to faster transport. 
This resembles noise-assisted transport found in many 
other systems&i^. The EST 7777(a), in contrast, increases 
for larger values of A, leading to slower transport to the 
drain. Thus, here the optimal transport efficiency is ob- 
tained in the purely coherent case (a = and A = 0). 
However, 7777(a) is always smaller than 777(a). We fur- 
ther observe that 7777(a) increases until a certain a c , after 



which it follows the curve of 777(a). The value of a c be- 
comes smaller with increasing dephasing rates. This hap- 
pens because the dephasing process destroys the coher- 
ences between nodes 1 and 2. Therefore after this point, 
the initial correlations do not significantly influence the 
transport properties anymore and 7777(a) w 777(a). When 
<p = 0, a tractable analytical expression for a c is possible: 



3 + 2A 



(2 + A3) 



(17) 



11 + 9A + 2A 2 

Our results show that initial correlations, induced by 
the source feeding the superposition state = (|1) + 
|2))/V% leads to faster transport than for feeding any 
other state of the form |^) = (|1) + e^|2»/>/2 with 
(j) ^ 0, even in the presence of dephasing. Additionally, 
if cos0 > one always has 7777(a) < 777(a), see Eqs. (fl3]) 
and ([HI). Therefore, when initial correlations are present, 
a smaller coupling to the environment (smaller values of 
a) is sufficient to avoid the dark state, compared to the 
situation without initial correlations. 



IV. SUMMARY 

In conclusion, we have shown that by connecting a 
source to the two end nodes of the V-shaped network it 
is possible to induce initial correlations between the cou- 
pled nodes and that these initial correlations can over- 
come the detrimental effects of the dark state, leading to 
complete transfer to the drain. The source can also feed 
the two end nodes independently, i.e., without initial cor- 
relations between the two nodes. Then, the dark state 
inhibits complete transfer. When increasing the coupling 
to the environment, the differences in the expected sur- 
vival times between the two types of sourcing processes 
diminish. We expect that the results obtained here also 
hold for larger networks that exhibit invariant subspaces, 
as for example described in^ for fully connected networks, 
inii for larger ring-like structures or inii for Erdos-Renyi 
graphs. Furthermore, our results are also related to the 
study of electrical currents through a network of two- level 
systems, since the current is related to the long-time limit 
of the time derivative of the EST. 
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